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ABSTRACT 


Preliminary experiments in the numerical prediction of large scale 
sea-Ssurface temperature anomalies are made using a 10-level primitive 
equation model with 300-km horizontal resolution covering a rectangular 
basin in the North Pacific. The model is first integrated over an 1ll- 
year period to statistical equilibrium using time dependent wind and 
thermal forcing. The monthly normal climatology, generated by the 
model in this long term integration, is then used along with observed 
sea-surface temperature anomalies to define the initial state. 

The results show that inclusion of salinity in the model may pro- 
duce more accurate predictions of high latitude cold anomalies. The 
most accurate of the predictions removed all of the anomalous diffusion 
terms in the governing equations. Data for the observed temperature 
anomalies were not available at depth, thus initialization of the 
anomalies below the surface was Superficial. Consequently, the effects 
of vertical advection were not fully realized, indicating a need for 
future experimentation with a more accurately defined initial profile 


of the anomalies in the initial conditions. 
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Ee. INTRODUCTION 


For today's meteorologist to make accurate, long range forecasts the 
understanding of the dynamics of atmospheric behavior is imperative. In 
the recent metecrological past, the importance of air-sea interactions 
has been realized and has given the forecaster another potential tool with 
which to improve his predictions. A better understanding of the inter- 
actions between the oceans and the atmosphere will also assist in ocean 
acoustic forecasting and other Naval operations. 

Like the atmosphere, the oceans experience seasonal, monthly, and 
daily variations in temperature and circulation patterns. Further investi- 
gation of these variations reveals the existence of anomalous behavior 
within the mean temperature and circulation patterns. It is believed 
that these variations influence directly the atmospheric weather patterns 
through the exchange of moisture, momentum, and heat across the air-sea 
interface. This thesis is concerned with the prediction of these 
anomalous oceanic temperatures and currents. 

During the past two decades it has become apparent through the studies 
of J. Namias and the late J. Bjerknes that the atmospheric response to 
sea-surface temperature anomalies may be worldwide. Bjerknes (1966, 1969, 
1972) has shown that large-scale, sea-surface temperature anomalies exist 
in the equatorial region and that these anomalies produce a response in 
the atmosphere which is not limited to that region alone. Likewise, 
Namias (1959, 1969, 1972b) discusses the possibility that such a response 
is also present in the middle latitudes. Furthermore, Namias (1969) 


Observed that during several years when intense sea-surface temperature 





anomalies were present in the North Pacific Ocean, there were distinct 
weather patterns over the United States. Because of the above pioneer- 
ing findings of Namias and Bjerknes, the North Pacific Experiment 
(NORPAX) was established in which a cooperative scientific effort is 
being made to understand the causes of these anomalies and their rela- 
tionship to the atmosphere. 

Numerical integrations of coupled ocean-atmosphere models have been 
reported by Manabe (1969), Bryan (1969), and Wetherald and Manabe (1972). 
By comparing the results of their coupled model with those of the ocean 
and atmosphere models integrated separately, these authors showed the 
Significant role ocean circulation plays in determining the state of 
the atmosphere. 

In turning our attention to the oceans, the first attempts at large 
scale ocean predictions was based upon the use of a prediction equation 
for the sea-surface temperature (Namias, 1965). In Namias' method, 
temperature T is a function of x,y,z, and time t. The local rate of 


change of temperature can be represented by 


OT _ oT Cee oe 
— = uae Vv ay W oe +g (1) 


where q is the rate of change of temperature from non-advective pro- 
cesses. Now if we let ( ) denote the long term normal state and (') the 


anomaly, then the monthly mean can be expressed by 


and 
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Then Eq. (1) takes the form: 


Ae ee eo (Tem!) = THT") 
oe ta") “x — (vtv") — 
- 0 (T+T! - 
- (wtw') oe einai o ke (2) 


Assuming the long term normal temperature satisfies the equation 


a - Va CHW +g (3) 


subtracting (3) from (2) gives the equation for the temperature anomaly: 


ones, 82, BT", x.) 8} 
oa -[{(utu') ax + (v+v') ay + (wtw') ve 
' OT ’ oT wor s 
-f[u eet By + ow ae sae (eee (4) 


Namias' first predictions (1965) retained only the advection of the 


oT 
—— + 


normal temperature by the anomalous currents, (u' ae eV, By)! algae (e 


(4). The anomalous currents (u', v') were obtained from the anomalous 


atmospheric wind stress observed during the period of the "forecast" 

and the long-term normal Ses ikes (T) from observations collected 

by the U. S. Weather Bureau. Since this prediction used data actually 
observed during the forecast period, it is not a prediction in the true 
sense and is referred to as a "Semiprediction". Although this simple 
method produced results which correlated positively with the observed 
temperature changes, Arthur (1966) demonstrated that additional terms 
should be included in the prediction equation. As a result, Jacob (1967) 
performed a semiprediction neglecting only the vertical advection terms 
and the anomalous heating. A Precis comparison of his results 


indicated an improvement in the prediction. 


i. 





More recently, Namias (1972a) has shown that the sea-surface tempera- 
ture patterns in the North Pacific can be predicted with some success 


using normal surface currents to advect the anomalies (neglecting all but 


= oc Be ona 
dy 


(Figs. 1, 2, 3, and 4), which will be discussed in detail later, show 


on the right side of Eq. (4)). The results of his study 


the importance of horizontal advection. These results plus the lack of 
temperature data at depth clearly indicates the need for a completely 
dynamical model. Therefore, it is the objective of this thesis to predict 
the evolution of the anomalous temperatures using a numerical model which 
contains all the advection terms in (4) and a parameterization of mixed 
layer dynamics as well. 

To make a comparative analysis, the period covered by’ this prediction 
will be the same as that used by Namias (1972a). Sea-surface temperature 
anomaly data were also obtained from Namias' 20-year normal observations. 
It should be noted, however, that this normal data is from the years 1947 
to 1966 and is not necessarily representative of the prediction period. 
Furthermore, due to the aforementioned lack of data at depth, we assume 
an anomalous temperature profile which decreases linearly from its value 
at the surface to zero at 150 meters. NORPAX scientists have shown that 
anomalies exist at depths greater than 150 meters; however, this depth 
was chosen to correspond more closely to Namias' study which assumed the 
temperature anomaly to be confined strictly to the surface. Further 
studies could be conducted varying the vertical profile of the tempera- 
ture anomaly to determine the importance of the initial condition. The 
initial anomalous current is geostrophic and its vertical average is 


zero. Unlike Namias‘ method which employed the use of observed 
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climatological data during the entire forecast period, this method uses 
model-generated climatology only to define the initial conditions. Con- 
sequently, following initialization the anomalous and mean fields are 


completely free to dynamically interact. 
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Pee tht (NUMERICAL WODEL 


A. BACKGROUND 

The basic model used in the following predictions is that of Haney 
(1974) as modified by Davies (1975). "The model is based upon the hydro- 
static and Boussinesq approximations. This means that the ocean is 
assumed incompressible, and the density is replaced by a conStant every- 
where except in the hydrostatic equation where it is multiplied by the 
acceleration of gravity. A major simplification is the neglect of 
salinity and the assumption that the density is a linear function of 
temperature atonen 


The equations which govern Haney's model are given by 




















du _ -1 ap u 2 97u 
dt —p acosd dA ey acosd> Oot a Deal pee’ ‘2 
oO OZ 
2 
me _-1 Op Ue . 2 dv 
dt pa dad a acoso a + aa Ca ad “) 
oO OZ 
gp 
i PG » (7) 
OW 1 du 3 _ 
ae acosd el dd eae os - 
2 
aT _ 2 oT 
a a Te 5 + 9 (7) ’ (9) 
Oz 
p = ee [Peewor = To) ] ’ (10) 


where the representation of the symbols is as follows: 


t time 
r longitude 
a) latitude 


il , 

Haney, R. L., "A Numerical Study of the Response of an Idealized 
Ocean to Large Scale Heat and Momentum Flux," J. of Physical Oceanog., 
Pe 4,°No. 2, p. 147, 1974. 
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height 

angular speed of the earth's rotation 
coriolis parameter 

radius of the earth 

acceleration of gravity 

eastward velocity component 

northward velocity component 

vertical velocity component 
temperature 

constant reference temperature 
density 

density of water at reference temperature 
coefficient of thermal expansion 


< fF Qn mDN 


= 


O 


ROUDUA HA 
O 


B. BOUNDARY CONDITIONS 

After Haney (1971) and Davies (1975), the determination of the upper 
boundary condition in the thermal energy equation is computed by calcu- 
lating separately the incoming solar radiation (Q.), the net upward long- 
wave radiation (Q,), the latent heat exchange (OQ), and the sensible heat 
exchange (Q.)- Together, these fluxes define the net downward heat flux 


as represented by 
= 0, aeewemon tO): (11) 


The back radiation SF 1s calculated from 


= 7 i772 7 2 4 
Q, (T) = -985(.39 - .05e " ) (1 - .6C°) oT. (12) 


where te is the sea surface temperature, = the atmospheric vapor 
pressure at 10 meters, C the fraction of cloud cover, and o the 
Stefan-Boltzmann constant. The exchange of latent heat and sensible heat 


across the ocean-atmosphere interface were computed from 


O., = =eezno Cc. lv] L fe it (13) 


E 


i> 





and 
= V Je = it 
om 0, OF |v| C, oo ae 


respectively, where 2a was the density of the air, oe the drag co- 
cifacient, |v | Pierce Specd at anemometer Jevel, L the latent heat 
of vaporization, a the atmospheric pressure, entre) the saturation 
vapor pressure at the sea-surface temperature, os the specific heat 

of air at constant pressure, and T. the atmospheric temperature at 


10 meters. Incoming solar radiation (Q,) was calculated with 


Q, = So(.74 - -6C) a) 


where oe was the solar insolation incident on a horizontal surface at 
the top of the atmosphere. 

The quantities So Pa oe Qe, tT. and |v| are prescribed functions 
of latitude (v depends on longitude--see below) and time taken from 
atmospheric climatological data. The fluxes Qn - ot , and Qu are 
additionally dependent on the computed sea-surface temperature, To a 
Therefore, an anomalous sea-surface temperature produces a corresponding 
anomalous heat flux which has the tendency to remove the temperature 
anomaly by virtue of the opposing sign in Eq. (11). 

"In addition to the thermal forcing, the ocean circulation was also 
forced by surface winds given on a time-dependent basis. Assuming an 
average geostrophic inflow angle of 10° and a frictional reduction of 
velocity by 10%, one could derive the surface stress from the surface 
geostrophic winds. If u and v_ are the eastward and northward surface 
winds, respectively, after inflow and frictional influences are considered, 
then 


Iv] = (a2 + v2)2/2 


16 





and the wind stress components are 


T 


- Pp. Cy |v| u 


T 
2 


2 
P. oe |v| a 
Therefore, the wind stresses are completely specified since u(y,t) 
and v(x,t) are completely independent of the ocean temperatures and 
currents in this experiment. Since the quantities which define the 
atmospheric forcing are derived from climatological] data, this thesis 
does not concern itself with the response of the ocean to anomalous 
atmospheric forcing. 

Lastly, the ocean bottom is assumed to be flat which requires the 
vertical velocity to vanish there. Thus integrating the continuity 
equation up to the surface and applying a free surface boundary condi- 
tion, one could obtain a prediction for the height of the surface. 
However, the present method assumes that the ocean surface 1s not a free 
surface but rather a balanced surface at which w = 0. Therefore, since 
w = 0 at the top and bottom, the vertical average of (8) shows that the 
vertical average current is non-divergent and consequently can be de- 
scribed by a streamfunction. The streamfunction is predicted from the 
vorticity equation which is derived from the vertical average of the 


equations of motion. 


C. RESOLUTION AND DATA 
TO improve the resolution, the 6-level model employed by Haney (1974) 


was expanded to include 10 levels (Davies, 1975). The vertically expanded 


Davies, R. W., A Numerical Parameterization of Wind Mixing ina 
Time-Dependent Baroclinic Oceanic General Circulation Model, M. S. Thesis, 
Naval Postgraduate School, Monterey, California, 1975, pp 56. 
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model now has fixed levels at 10, 30, 60, 100, 150, 225, 350, 700, 1500, 
and 3000 meters. 

As mentioned previously, to compare the results with those of Namias 
(1972a), sea-surface temperature anomaly data covering the same period 
were obtained from J. Namias of Scripps Institution of Oceanography. 

The surface temperature data thus obtained covered only the area 20°N to 
60°N, 130°E to 110°W in the North Pacific and the grid spacing intervals 
for the latitude and longitude were 5° and 5° respectively. Therefore, 
Namias' data were expanded to the equator and to 65°N by extrapolating 
the anomalous data to zero at the most distant boundaries. Additionally, 
the data over the entire domain were interpolated to fit the 2.0° lati- 
tude and the 2.8° longitude grid spacing of the present model. 

As a result of the expansion, the predictions outside the original 
data field should be ignored. Consequently, for the convenience of the 
reader, an overlay (Fig. 5) with the locations of ocean vessel stations 
November (N), Papa (P), and Victor (V) as well as the site of the NORPAX 
"Pole" experiment, has been superimposed on the normal fields and each 
related prediction. The superficial land mass boundaries of the NORPAX 
grid are represented by the shaded area and the original data field by 


the area between the dashed lines. 


D. INITIAL CONDITIONS 

To completely define the initial state of the ocean for dynamical pre- 
dictions, a considerable amount of model climatology was used along with 
the available observed data. The model climatology was generated by in- 
tegrating the model over an ll-year period starting from a stratified 


initial state and using atmospheric data which contained the observed 
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normal seasonal data. This model climatology was used for initializa- 
tion and verification as described below. To provide continuity in 

the atmosphere wind and thermal forcing, the annual mean and amplitude 
and phase of the first two harmonics were calculated from atmospheric 
climatological data as defined by Davies (1975). The atmospheric data 
were averaged along each grid latitudinal line, except for v, the north- 
ward component of the surface geostrophic wind, which was averaged 

along each meridian. 

The resulting initial normal currents and the 30, 60, 90 day normal 
currents at level 2 (30 m) are shown in Figs. 6, 8, 10, and 12, where 
"initial" is defined as October 15th of year 11. As seen in Figs. 6 and 
12, there is an almost negligible change in the normal ae during 
that period. The corresponding initial, 30, 60, and 90 day normal tem- 
perature fields are shown in Figs. 7, 9, 11, and 13. Although there is 
little change in the temperature gradient in the middle latitudes, the 
southward movement of the main baroclinic zone indicates a cooling trend 
which was expected during the winter months. 

To begin the predictions, the total field must be defined in the 
initial conditions. Thus, the normal part of the initial conditions is 
obtained from the model-generated climatology while the anomalous part 
1s obtained from observations. The observed sea-surface temperature 
anomalies are extrapolated linearly from observed values at the surface 
tO zero at Tevel 5 (150 meters). This rather arbitrary procedure of 
defining the anomalous temperatures below the surface is done because 


of the lack of sub-surface temperature data. 
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The total current (normal plus anomaly) can be represented by 


ra * 


ee toa 7) 


O 
* 
where () is the vertical average current, oa an Vaz, )\1as-the 
“A ‘ ~H 
mertical shear current, [( )-( )], and H is 


the total ocean depth. Expanding (17) into its normal and anomalous 


parts one obtains 


A an 


* - 
W=wWwtw =wtwtwey' (18) 


~ * 
where \V and \V are the normal part of the vertical average current 


and the normal part of the vertical shear current respectively. The 
anomalous part of the vertical average current WW") is assumed to be 
zero and the anomalous part of the vertical shear current av") is the 
geostrophic current calculated from the above temperature anomaly using 
the thermal wind relation. Since 


g PF 


taking the partial derivative with respect to z produces 


I 


a\V 1 
= Df Ik x V (-9q) ° 
O 
Substituting (10) for 0 
P ape Ger + constant , 


where the gradient of PP becomes 


VO Soe rg NMP : 
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Therefore, 


a Sik x VT (thermal wind relation). 


It is now assumed that the anomalous part of the vertical shear current 


satisfies this relationship, 


) 
Pr. Se is x Vane. 


Then integrating from -H to some depth 2z 


* * ag e 
W'(z) = \v'(-H) + rz 1k x vf Peace (19) 
=r 
* * 
where \V'(-h) is unknown. If the vertical average of \V' is taken, 
then this quantity can be set equal to zero by definition. Thus, from 


(9) 


Oo 


A oO 2 

* 2 f \ ag 1 

ao = ey aa asks rim t = 

Wea Wwds + Bex ve f cf raryazso, 
-H 


=) 9 =H 


Or 


I 
o>) 
e 


©) = . 
W'(-H) + = Ik x vif f Pacyaz = (20) 
-H — 


* 
Then (20) can be solved for \y'(-H) and substituting into (19) to give 


the result 


bs Oo Lb, 
ae ag : i 
V weBaevi f e- 5 [ ( f) rrac)az) : (21) 
-H —H =F 
a el 
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In Eq. (21) f£ is not allowed to become smaller than its value at 
13.25°N. This geostrophic approximation to the vertical shear current 
anomaly is then used as the initial condition in the present model. 
Having defined the initial temperature and currents by combining 
the model generated normal fields with the observed anomaly fields, the 
model is run for 90 days of simulated time. During this time, the total 
temperatures and currents interact and evolve together. At any given 
time later (e.g., 60 days) the predicted anomaly is defined to be the 
total field predicted for that time minus the model generated normal 


for that same time. 


Ee WHE PREDICTION EQUATION 

Since several cases were studied uSing variations of the thermodyna- 
Mic prediction Eq. (9) which governs the total temperature field, a 
detailed analysis of that equation follows. 

Neglecting the convective adjustment (0 (T)) in the analysis, (9) 


can be written 





2 
Gi on 2 ar 
a, ae Selina (22) 
OZ 
Expanding (22) 
2 2 2 
T fy At 
Stu tvet twee a ac fey Sa OS (23) 
Y d dy Oz 
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and substituting further for the normal and anomalous parts, (23) takes 


the form 
am oT’ ..—-oT - oT’. , Or ecm noabor <= oT" 
Ze a = -[{ ne fF ua eu ae ae an + SS +: Ye 
ee ee, ge Se yt Shy ye SEY 
oy y 2 Zz 3 é 
2- ae. 2- Za 2— 2m. 
sapet+ oo ee (24) 
Ox Ox oy dy Oz Oz 


Setting T', u', v', and w' equal to zero, then the normal field can be 


represented by 


me (Gee oe? ta ta: (25) 


Subtracting (25) from (24), the equation governing the anomalous tempera- 


tures becomes, 





cae. ee age, 5 OT’, BT, Ee ce 
ae Cae Ox 3 , dy , Oz 
oa, a, Zee, 
Oz OZ Ox Oy Oz 
The following investigation will reveal that initially \y'*VT' = 0 


by analysis of (21). If we let "D" represent the depth at which the tem- 
perature anomaly becomes zero as defined earlier, then the temperature 


anomaly (T') at some level z< 0 is 


Tamer vO) me deen) es TK y+ 0} le 
D -Z 
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where T'(x,y,0) is the sea-surface temperature anomaly and T'(x,y,-D) 


is the temperature anomaly at D = 150 meters. Since T'(x,y,-D) 


definition, then 
Zz 
fhe — lb (x,y, 9) cL a = * 


Substituting (27) into term (a) of (21) gives 


Zz 
I fees iz) © (ny, 0) ; 
-H 
where 
Z 
aie) = f a + 2yaz 


i 


and (tf) is a dummy variable representing the vertical coordinate. 


wise, (28) can be substituted into term (b) of (21) to produce 


\e) zZ Oo 
= f [ diggol eb Al Buz) Py , 9) dz 


=H = =i 


feel (x, 0). 
Therefore, (21) becomes 


* “Aw 
Vie = ik x Vi (F(z) T!(x,y,0) - FT (x,y,0)] 


Br by tactowang T"' (x,y,0) 


‘i See 2) Vaiaetyy/ 0) (F(z) - Ft . 
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(27) 


(28) 


Like- 


(29) 





“Aw 
From this, it is clear that F(z) and F are not functions of x and y and 


can be brought outside the "del" operator allowing 


+ ag a ' 
v= FZ «(CF (2) - F) lk x aes 27 5:0) -< (30) 


Consequently, 


x 


ies » _ OF _P ' ; i . 
\V ae - (F(Z) F) (tk x vee cary 0) Vat (x,y,0) (1 + a 


_ a _P Be ' ; ' 
= —p (F (2) F) (1+ D) [Ik pr Vt Cy, 0} ei Cy 0) 1 (315 


ee atl 
a b 


and since term (a) of (31) is parallel to the anomalous temperature 
field and term (b) of (31) is perpendicular to the same field, then 


* wn 
ew T' = 0. Since \y' = 0, we have 


a SO 


which allows wre and ve to be neglected in (26). Figures 14 

and 15 clearly show that initially the anomalous currents are geostrophic 
and parallel to the anomalous isotherms. It should be noted that these 
currents need not remain parallel to the anomalous isotherms or geos- 
trophic during the prediction period; however, careful examination of 

the predictions indicates they do. It is not surprising that this type 
of motion remains geostrophic, but it 1S surprising that it remains equi- 


valent barotrophic. Therefore, for the purpose of this interpretation, 


we will assume that \y' ° vee is negligible during the forecast period. 
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Additionally, 


i] 
wo can be further evaluated. From the continuity 
equation (8), 
Seu WHO, 
then 
ow — ow 
—— a + — + e 1 
5 + bre WV 7 Se \V 0 
where 
oe Me -W=0 F 
implying 
one 
OZ 


+ ah) = 
ce \V 0 
Thus integrating up from 


the level z to the surface z 


— O ; 
oO 
w' = ¥ Va e \W' dz + w' (0) (32) 
z 
where w'(0) = 0 as defined by the boundary conditions. Therefore, (32) 
can be written 
O 
vo f ie av as. , (33) 
Z 
“A * 
and since Me av = OC, es ei = oe -\' allowing the substitution of 
(9@) into (33) as 
p F) 
ie Fey —F 
Ww ag f Mi [ 
Zz 


(34) 
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Samce thercorrtolis parameter (f) is a function of y, and T'(x,y,0) is 


not a function of z, Eq. (34) becomes 


ee » (eg ' = 
wi = Me le x se PCO) EZ) ] (35) 


where 


O 
F(z) -f (F(z) - F) dz. 
z 


Hence, evaluating (35) 


gee [Vge(— SS SE Ye OD 5 ue Og ST V1) 5) F(z) ] ; (36) 


If we denote the geostrophic current, integrated from the level z to the 


surface z = 0 by (Ug, Vg) , 


Ug (z) => = BS yO) F(z) 


and 


wee) = OG OR 8). (2) : 


then (36) can be written 


re) Pe) 
wi = [sx (Ug) + ae ED : (373 
Therefore, 
eee oes ee ixey, 0) =, _ _ 8 
o> 2 oy ax PO) ll re CA 
and so, 
oT! B oT! 
| elie Ia eh il 
wl = Vg (z) a (38) 
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The approximation for the vertical velocity w' applies after integration 
begins as well as initially. This result shows that w' is the vertical 
motion produced by the divergence of the goestrophic current. Thus, a 
northward moving layer of surface water (Vg > 0) converges to produce 
sinking motion below the surface (w' < 0), while southward moving surface 
water (Vg < 0) diverges and produces rising motion below the surface 

(w' > 0). As a result of the above evaluated terms of (26), the final 


prediction equation for the anomalous temperatures becomes 


1 2 3 4 5 6 

ieee, Sf ye gee ee moles Ol 6 oT 
Pra is Ox va Oy xO ae BS £ Vg (2) Oz 

t 2nt 2A 2 

SG) |) Ae ee (39) 
fe dz 2 2 2 
Ox dy Oz 
7 8 9 


The first two terms on the right-hand side of Eq. (39) represent 
the horizontal advection of the anomalous temperature by the normal 
current. These are the only terms used by Namias (1972a) to make his 
prediction. Terms 3 and 4 are ane horizontal advection of the normal 
temperature by the anomalous current. Term 5 is the vertical advection 
of the anomalous temperature by the normal current. The vertical advec- 
tion of the normal temperature by the anomalous current and the vertical 
advection of the anomalous temperature by the anomalous current are 
terms 6 and 7 respectively. Term 8 is the horizontal eddy diffusion of 
the anomalous temperature and term 9 is the vertical eddy diffusion of 
the anomalous temperature which also includes the anomalous part of the 


surface heat flux (see discussion below (15)). 
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The motivation behind the inclusion of all the terms of (39) is, of 
course, to make a more accurate prediction. As a result of Jacob's (1967) 
work, which showed an improvement over Namias' predictions (1972a) by 
including terms 3 and 4 of (39), and since vertical advection is so im- 
portant in the theory of the permanent thermocline, the implication is 
that all the advection terms should be included when making a prediction 


of the sea~surface temperature anomaly. 
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ITI. RESULTS 


Before discussing the results of the predictions generated by the 
complete dynamical model, a brief description of the observed anomalous 
temperature fields and the results of Namias'(1972a) predictions fol- 
lows. The initial, 30, 60, and 90 day observed anomalous temperature 
fields, as seen in Figs. 15, 16, 17 and 18, are averages over three 
months centered about the months of October, November, December, and 
January respectively. That is, the initial anomalous temperature field 
is the average of the anomalous temperatures for the months September, 
October, and November (i.e., the Fall season) of 1971. 

From Figs. 15 through 18 two cold anomalies and one warm anomaly 
appear to dominate the area of interest. The cold anomaly near the 
western boundary is observed to decrease in amplitude almost linearly 
during the first 60 days with a small southward movement. The double 
cell cold anomaly near the northeast boundary decreases in amplitude 
in the western part but maintains its relative amplitude and position in 
the eastern part. Finally, the warm anomaly near the center of Fig. 15, 
increases in amplitude with the leading edge of the maximum anomalous 
temperature contour moving eastward some 13° longitude. 

If one observes Namias' 90 day prediction of the warm anomaly (Fig. 
4), aS described above, the easterly movement of the anomaly, as 
advected by the normal winter sea-surface currents (Fig. 3), is a good 
approximation. However, the amplitude of the anomaly in his prediction 


remains relatively constant, when it was obServed to increase. In an 
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Sttempet tO explain Chis discrepancy, Namias calculated the surface heat 
fluxes, but the results obtained indicated a negligible (and at times 
opposed) role of these fluxes. It is therefore suggested that the over- 
simplified use of the prediction equation, as described earlier, and the 
lack of data at depth is responsible for not predicting accurately the 
amplitude of the anomaly. 

Attempting to improve the predictions made by Namias (1972a), the 
completely dynamical model of this thesis was run for a corresponding 
90-day period using the initial anomalous temperature field (Fig. 15) 
and the model generated climatological data along with the boundary con- 
ditions as described earlier. Four cases uSing variations in the 
boundary conditions and this prediction equation are discussed below. 

A brief description of the four cases will be given first, followed 


by a qualitative discussion of the results obtained. 


GASE 1: 
In CASE 1 the anomalous temperature prediction contained all the 
terms implied by Eq. (39). The results of CASE 1 are shown in Figs. 19- 


24. 


CASE 2: 

In CASE 2 the model calculated the 90 day prediction as in CASE l 
except in the surface thermal boundary condition, Eqs. (12), (13), and 
(14), the normal climatological sea-surface temperatures were substituted 
for the predicted sea-surface temperatures. This effectively removed the 
anomalous part of the surface heat flux. The results of CASE 2 are shown 


imeeas. 25o=30. 
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GASE 3: 
CASE 3 is similar to CASE 2 except normal climatological temperatures 
were substituted in place of the predicted temperatures in the diffusion 


terms of (9). This effectively eliminates the anomalous diffusion terms, 
—- Dit ae 

facet — 2 —_ ]) , from (39). The results of CASE 3 are shown in 

ox Oy OZ 


Migs. 31-36. 


CASE 4; 
CASE 4 is similar to CASE 3 except normal climatological currents 


(u,v) were substituted for the predicted currents in the diffusion terms 
2 2 
Ce 

2 2 
OZ OZ 5 
This has the effect of removing the anomalous diffusion terms (A V y', 

ee Aaa 
K emu and A Vev', K apy 
p 2 

OZ OZ 


currents. The results of CASE 4 are shown in Figs. 37-42. 





(A Bio yk Mor thewpGimitive Eqs. <5) and.(6).. 





) from the equations governing the anomalous 


Concentrating on the results of the anomalous temperature predictions 
of CASE 1 (Figs. 15, 20, 22, and 24), one observes that all the anomalies 
appear to decrease in amplitude with the cold anomalies disappearing com- 
pletely by the end of that 90 day prediction. The amplitude of the warm 
anomaly initially was found to be .6°C located at 168°W, 40°N and is 
reduced to .1°C after 90 days with its center at 162°W, 37°N. Asa 
result of the above observations, cases 2, 3, and 4 were designed to 
determine the cause of the amplitude decrease or disappearance of these 
anomalies. Furthermore, the western edge of the warm anomaly is observed 
to move eastward while the eastern edge remains relatively stationary. 
Below 20°N some anomalous activity appears to have been generated; 
however, since this lies outside the area of initial observed data, it 


is disregarded and further discussion of this region will be omitted. 
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As in CASE 1, the cold anomaly at the northeast boundary disappears 
Pieensh 2d me 20,020, and 30). Analyzing the initial vertical 
profile of the total temperature field, it was found that the lapse rate 
was unstable thereby implying that convective adjustment (term 6 (7) sya 
(9)), is the cause of the cold anomaly vanishing. This rapid disappear- 
ance of the cold anomaly immediately raises the question of the impor- 
tance of salinity. If salinity were included in the present model, it 
is speculated that the cold anomaly could maintain itself. That is, the 
cold anomaly at the surface would not necessarily be unstable if a (nega- 
tive) salinity anomaly were present to stabilize the water in the column. 
Another possibility is that salinity is negligible, but the cold anomaly 
extends all the way to the ocean floor. However, this explanation seems 
less likely since it indicates the presence of such a large anomaly of 
heat content. 

Investigation of the initial vertical temperature profile of the 
total field in the vicinity of the cold anomaly near the western boundary 
reveals that it too is convectively unstable. However, this anomaly damps 
less rapidly in CASE 2 than does the cold anomaly at the northeast 
boundary which would indicate its decrease in amplitude is not entirely 
due to convective adjustment. Additionally, the amplitude of the central 
warm anomaly was found to decrease less (to a value of .2°C) than in 
CASE 1. Except for the above two differences, the 60 and 90 day predic- 
tions of cases 1 and 2 are relatively similar. Therefore, with respect 
to the observations of the first two cases, it is concluded that the 
anomalous surface heat flux induced by the anomaly itself is not the 
major cause of the damping of the anomalies, and that convective adjust- 
ment is the major cause of the cold anomalies vanishing. The effects of 


diffusion will be discussed in cases 3 and 4 below. 
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Figures 15, 32, 34, and 36 show the results of CASE 3, in which there 
is no eddy diffusion of the temperature anomaly. It is seen that the 
cold anomaly at the northeast boundary still disappears as in cases 1 and 
a, the een cold anomaly, however, appears to maintain its magnitude 
for a greater duration than in the previous two cases; the greatest dif- 
ference occurring during the 60 to 90 day prediction period. The western 
edge of the warm anomaly moves eastward about the same distance as before 
but its shape is aiieered and the eastern edge remains unchanged as in 
cases 1 and 2. The major difference observed in CASE 3 is that the maxi- 
mum amplitude center of the warm anomaly only decreases from .6°C ini- 
tially to .4°C at the end of the 90 day prediction. Consequently, the 
results of CASE 3 imply that the anomalous diffusion terms of (39) were 
Significant in producing the heavy damping in cases 1 and 2. Since re- 
moving the anomalous diffusion terms of (39) improved the results of the 
prediction considerably, but did not completely eliminate all of the 
damping, CASE 4 was devised in which the diffusion of the anomalous cur- 
rents was also removed. The removal of these additional diffusion terms 
should eliminate some damping because of the dependence of the tempera- 
ture on the currents through the geostrophic relationship as discussed 
previously. 

With respect to dissipation, movement, and amplitude of the anomalies, 
CASE 4 (Figs. 15, 38, 40, and 42) produced results nearly identical to 
those of CASE 3. Therefore, it 1S concluded that the diffusion terms 
governing anomalous currents in the primitive equations are seemingly 
unimportant in the prediction of sea-surface temperature anomalies. Never- 
theless, CASE 4 represents a more consistent model in that there is no 


diffusion of either the anomalous currents or temperature. 
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If one now compares the results of Namias' prediction (Fig. 4) with 
the results found in this study using a completely dynamical model, it 
is first observed that the 0° isotherm of the warm anomaly has about the 
same configuration and location. Secondly, the location of the center 
of maximum amplitude in the predicted warm anomaly is much closer to the 
observed 90 day position in Namias' study than in this study. However, 
if one closely examines the location of this anomaly center in Figs. 15 
and 16, it is seen that the anomaly center moves eastward almost 15° of 
longitude during the first 30 days and then remains relatively station- 
ary. This observation would imply that the movement of the anomaly's 
center is not due simply to advection alone as indicated in Namias' pre- 
diction. Lastly, the amplitude of the predicted warm anomaly decreases 
in all cases of this study but the amplitude of the same anomaly in 
Namias' prediction remains constant as expected when employing horizontal 
advection alone. Therefore, the vertical advection in this model plays 
an important role in the evolution of such an anomaly. Whether that 
role in the present prediction or any future predictions is simulated 
eeeuralely by the model depends critically on the vertical profile 


assumed for the anomaly. 
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IV. CONCLUSIONS 


This study has shown that the method of prediction usSing model gener- 
ated climatology along with observed temperature anomalies ina fully 
dynamical numerical model of the ocean circulation can be successful. 
Also, the simple geostrophic initial conditions for the anomalous cur- 
rents is successful away from the equator. However, the quality of the 
above results depends on the precise initial data available. 

In view of this study, the need for temperature data (normal and 
anomalous)in the mixed layers of the ocean is concluded to be the dominat- 
ing factor influencing the prediction of the sea-surface temperature 
anomalies. This appears quite evident from our observation discussed 
above in reference to Figs. 15 and 16. Since this rapid movement of 
the anomaly center cannot be explained by horizontal advection alone, 
1t 1S speculated that a sub-surface warm anomaly may have existed down- 
stream which was brought to the surface. Thus, vertical advection may 
have played a significant role; however, the lack of data beneath the 
surface makes these speculations virtually impossible to investigate. 

The present model, with no salinity, cannot accurately predict the 
evolution of high latitude cold anomalies which do not extend to the 
ocean floor. Therefore, initialization should include salinity; however, 
at present, large scale anomalous salinity data is not available. 

When comparing the results of cases 3 and 4, it is observed that 
the removal of the diffusion terms governing the anomalous currents had 


Pvttle etfect, if any, on the 90 day prediction. Thus, further 
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predictions using a dynamical model should, for consistency, omit all of 
the anomalous diffusion terms (CASE 4). 

Although the results of this study did not produce predictions which 
corresponded exactly with the observed anomaly field, the use of a fully 
dynamical model in making these predictions has given us a greater in- 
Sight as to the role played by vertical advection and diffusion. It is 
therefore believed that further studies should be conducted using this 
fully dynamical model with more accurate initial data and with salinity 
included. 

Additionally, this study uses only climatological data to define the 
atmospheric forcing. However, an anomalous atmosphere would produce 
anomalous heating and wind stresses which could produce anomalous features 
in the ocean. For example, Namias (1965) investigated the effects of an 
anomalous wind stress and showed that this anomalous atmospheric forcing 
may be important in some cases. Thus, further investigations should be 
conducted to determine the response of the ocean to various anomalous 


atmospheric forcing. 
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